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ABSTRACT 



Our universe is observed to be accelerating due to the dominant dark en- 
ergy with negative pressure. The dark energy equation of state (w) holds a key 
to understanding the ultimate fate of the universe. The cosmic voids behave 
like bubbles in the universe so that their shapes must be quite sensitive to the 
background cosmology. Assuming a flat universe and using the priors on the 
matter density parameter (Q m ) and the dimensionless Hubble parameter (h), 
we demonstrate analytically that the ellipticity evolution of cosmic voids may 
be a sensitive probe of the dark energy equation of state. We also discuss the 
parameter degeneracy between w and Q m . 

Subject headings: cosmology: theory — large-scale structure of universe 



R ecent observation s have revealed that ou r universe is flat and in a phase of acceler- 
ation (IRiess et al.l 1 19981 ; iPerlmutter et al.l Il999l ; ISpergel et al.l 120031 ). It implies that some 
mysterious dark energy fills dominantly the universe at present epoch, exerting anti-gravity. 
The nature of this mysterious dark energy which holds a key to understanding the ultimate 
fate of the universe is often specified by its equation of state, i.e., the ratio of its pressure to 
density: w = Pde/ Pde- The anti-gravity of the dark energy corresponds to the negative value 
of w. The simplest candidat e for the dark energy is the vacuum energy (A) with w = —1 
that is constant at all times (lEinsteinlll917l). Although a l l current data are consistent with 



the vacuum energ y model (e.g., IWang fc Tegmarkl 120041 ; iJassal et al.l 120041 ; iPercivall 12005 



Guzzo et al. 



2008) , the notori ous failure of the theoretical estimate of the vacuum energy 



density (see lCaroll et al.lll992l . for a review) has led a dynamic dark energy model to emerge 
as an alternative. In this dynamic dark energy models which is often collectively called 
quintessence, the dark energy is described as a slowly rolling scala r field with time- varying 
equation of state in the range of — 1 < w < (jCaldwell et al.lll998l ). 
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The following observables have so far been suggested to dis criminate the da r k en- 
ergy models: the l uminosity-distance m easure of type la supernova (IRiess et al.l 12004 120071 ; 



Davis et al 



2007 



Kowals 



ri et a. 



mass (IWang &: Steinhardt 



199 



200 81); the abundance of galaxy clus ters as a function of 



Haiman et al 



tic osci 



Coorav 


2004 


Huterer 


2001 



lation s in the galaxy power sp ectrum (IBlake fc Glazebrook 



2001 



Weller et al. 



20021), the baryonic 



acous- 



2003; Hu fc Haiman 



Takada fc Jain 



2004 



Seo fc Eisenteinl 120051). and the weak g ravitational lensing effect (IHu 



2003 



1999 



Song fc Knox! 120041 ) . True as it is that these observables 



can constrain powerfully the value of w, it is still quite necessary and important to find out 
as many different observables as possible for consistency tests. 

Another possible observable as a dark energy constraint may be the shapes of the cos- 
mic voids. As the voids behave like bubbles due to their extremely low densities, their 
shapes determined by the spatial distribution of the void galaxies tend to change sensitively 
according to the competition between the tidal distortion and the gravitational rarefaction 
effect. Therefore, the shape evolution of the voids must depend sensitively on the background 
cosmology. In this Letter we study the ellipticity evolution of cosmic voids in the QCDM 
(q uintessence + cold dark matter) model with the help of the analytic formalism developed 
by iPark fc Led (120071 ) and explore the possibility of using it as a complimentary probe of 
the dark energy equation of state. 



According to IPark fc Led (120071 ) , the shape of a void region is related to the eigenvalues 
of the local tidal shear tensor as 

1 + (5 V - 2)v 2 + /i 2 
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where {Aj}f =1 (with Ai > A2 > A3) are the three eigenvalues of the local tidal field smoothed 
on void scale, S v is the density contrast threshold for the formation of a void: S v = Y^=i 
and {/!, u} (with v < fi) represents a set of the two parameters that quantify the anisotropic 
distribution of the void galaxies. They defined the void ellipticity as e = 1 — v and evaluated 
its probability density distribution as 

p(l-e;z) = p(v;z,R L )= p\p,,u\6 = S v ;a{z,R L )]dfj, 
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Here, cr(z, Rl)) = D 2 ( z ) A 2 {k)W 2 {kR Relink is the linear rms fluctuation of the matter 
density field smoothed on a Lagrangian void scale of Rl at redshift z where D(z) is the 
linear growth factor, W{kRi) is a top-hat window function, and A 2 (k) is the dimensionless 
linear power spectrum. Throughout this study, we adopt the linear power s pectrum of the 



cold dark matter cosmology (CDM) that does not depend explicitly on w (IBardeen et al. 
1986h . 



Equation ([3]) was originally derived under the assumption of a ACDM model (w = — 1). 
We propose here that it also holds good for the case of a QCDM (quintessence+CDM) model 
where the dark energy equation of state changes with time as w(z) = wq + w a z/(l + z) 



(jChevallier &: Polarskill200ll ; lLinderll2003l ) where Wq is the value of w at present epoch and 



w a quantifies how the dark energy equation of state changes with time. Then, we employ 



( Basilakos 


2003; 


Percival 


2005) 
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The CDM density parameter Q m and the dark energy density parameter Qq evolve with z 
respectively as 

Q m0 {l + z) 3 



O m (^) 



QO 



(9) 



E 2 {z) ' " vv ~ y E 2 (z)(l + z)f( z )' 
where Q m o and Qqo represent the present values. Equation (jBJ) implies that the mean elliptic- 
ity of voids decreases with z. A key question is how the rate of the decrease changes with the 
dark energy equation of state. Since most of the recent observations indicate t hat the dark 



energ y equation of state at present epoch is consistent with w = —1 (e.g., see iGuzzo et al. 
20081 . and references therein) we focus on how the mean void ellipticity depends on the value 
of w a . Even in case that Wo = — 1, if w a is found to deviate from zero, it would imply the 
dynamic dark energy, disproving the simple ACDM model. 

To explore how the void ellipticity evolution depends on w a , we evaluate the mean 
ellipticity of voids as e(z) = e p(e; Rl, z)de for different values of w a through equations 
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©- (JS]). The other key cosmological parameters are set at fl m = 0.75, Qq = 0.75, h = 0.73, 
cr 8 = 0.9 and u> = — 1. When the abundance of evolution of galaxy clusters is used to 
constrain the dark energy equation of state, the cluster mass is usu ally set at a certain 



thres hold, Mr, defined as the mass within a certain comoving radius (I Wang fc Steinhardt 



19981 ). Likewise, we set the Lagrangian scale of a void, Rl at 4/i _1 Mpc, which is related 
to the mean effective radius of a void as Re = (1 + 5 v ) -1 ' 3 .Ri,/(l + z). The Lagrangian 
scale Rl = 4/i _1 Mpc corresponds to the mean effective size of a void at present epoch, 
R E ~ 8.5/i- x Mpc. 

Figure [1] plots e(z) for the four different cases: w a = —1/3, 0, 1/3 and 2/3 (long-dashed, 
solid, dashed, and dotted line, respectively). As can be seen, the higher the value of w a is, 
the more rapidly e(z) decreases. It also suggests that e(z) is well approximated as a linear 
function of z in recent epochs (0 < z < 0.2). Therefore, we fit e(z) to a straight line as 
e(z) ~ A v z + B v . Varying the value of w a in the range of [0, 2/3], we compute the best-fit 
slope A v . The range, < w a < 2/3, corresponds to the dark energy equation of state range, 
— 1 < w < —0.9. The result is plotted in Fig. [2j As can be seen, the void ellipticity evolves 
more rapidly as the value of w a increases. That is, the void ellipticity undergoes a stronger 
evolution when the anti-gravitational effect is less strong in recent epochs. Note that A v 
shows a noticeable 30% difference as the dark energy equation of state changes w from —1 
to -0.9. 

We have so far neglected the parameter degeneracy between w and the other key pa- 
rameters. However, as the dependence of the void ellipticity distribution on the dark energy 
equation of state comes from its dependence on A 2 (k; Q m o, c% h, w), it is naturally expected 
that there should be a strong parameter degeneracy. Here, we focus on the degeneracy be- 
tween Q m0 and w. First, we recompute A v , varying the values of Q m0 and w with setting 
w a = 1/3. The left panel of Fig. plots a family of the degeneracy curves in the fl m0 -Wo 
plane for the three different values of A v . As can be seen, there is a strong degeneracy 
between the two parameters. For a given value of A v , the value of wq increases as the value 
of Q m o decreases. A similar trend is also found in the fl m o-w a degeneracy curves that are 
plotted in the right panel of Fig. [3] for which the value of Wq is set at —1. It is worth noting 
that this degeneracy trend is orthogonal to that found from the cluster abundance evolution 
(see Fig. 3 in Wang & Steinhardt 1998). Thus, when combined with the cluster analysis, the 
void ellipticity analysis may be useful to break the degeneracy between Q m o and w. 

We have shown that the void ellipticity evolution is in principle a useful constraint of the 
dark energy equation of state. We have also shown that it provides a new degeneracy curve 
for Q m0 and w. When combined with the cluster abundance analysis, it should be useful 
to break the degeneracy. Furthermore, unlike the mass measurement of high-;? clusters 



- 5 - 



which suffers from considerable scatters, the void ellipticities are readily measured from the 
positions of the void galaxies without requiring any additional information. 

To use our analytic tool in practice to constrain the dark energy equation of state, 
however, it will require to account for the redshift distortion effect since the positions of the 
void galaxies are measured in redshift space. In our companion paper (Park fc Lee 2009 in 



prepa ration), we have analyzed the Millennium Run Redshift-Space catalog (jSpringel et al. 



20051 ) and determined the ellipticity distribution of the galaxy voids. From this analysis, it 
is somewhat unexpectedly found that the void ellipticity distribution measured in redshift 
space is hardly changed from the one in real space. In fact, this result is consistent with 
the recent claims of Hoyle & Vogeley (2007) and that of van de Weygaert (2008, private 
communication) who have already pointed out that the redshift distortion effect has only 
negligible, if any, effect on the shapes of voids. We hope to constrain the dark energy 
equation of state by applying our theoretical tool to real observational data and report the 
result elsewhere in the near future. 
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Fig. 1. — Mean ellipticity of the voids with R L = Ah 1 Mpc as a function of z. 
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Fig. 3. — Contours of A v in the f2 m0 -u>o (left) and in the Q m -w a (right) plane. 



